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Accurate results are obtained for the low temperature magnetization versus magnetic field of 
Heisenberg spin rings consisting of an even number N of intrinsic spins s — 1/2, 1, 3/2, 2, 5/2, 3, 7/2 
with nearest-neighbor antiferromagnetic (AF) exchange by employing a numerically exact quantum 
Monte Carlo method. A straightforward analysis of this data, in particular the values of the level- 
crossing fields, provides accurate results for the lowest energy eigenvalue En{S, s) for each value of 
the total spin quantum number S. In particular, the results are substantially more accurate than 
those provided by the rotational band approximation. For s < 5/2, data are presented for all even 
N < 20, which are particularly relevant for experiments on finite magnetic rings. Furthermore, we 
find that for s > 3/2 the dependence of En{S,s) on s can be described by a scaling relation, and 
this relation is shown to hold well for ring sizes up to N = 80 for all intrinsic spins in the range 
3/2 < s < 7/2. Considering ring sizes in the interval 8 < N < 50, we find that the energy gap 
between the ground state and the first excited state approaches zero proportional to 1/N a , where 
ct ~ 0.76 for s = 3/2 and a w 0.84 for s = 5/2. Finally, we demonstrate the usefulness of our 
present results for En{S,s) by examining the Fei2 ring-type magnetic molecule, leading to a new, 
more accurate estimate of the exchange constant for this system than has been obtained heretofore. 

PACS numbers: 75.10.Jm, 75.50.Ee, 75.40.Mg, 75.50.Xx 



I. INTRODUCTION 
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Since the early 1990s, the field of magnetic molecules 
has blossomed, and the number of different species 
that exist is increasing rapidlyi^iiii^ In particular, 
there is a large family of so-called ring-type magnetic 
moleculea 3 ! 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 : 16 that we focus on in the 
present work. Within such molecules there are embed- 
ded an even number N of identical paramagnetic ions of 
intrinsic spin s occupying N equally-spaced sites defining 
a ring. Each such ion ("spin") is coupled to its two near- 
est neighbors via an AF exchange interaction, resulting 
in systems that can often&LifiiiLiSiM be well represented 
by an isotropic Heisenberg model with a single exchange 
energy, J > 0, of the form 
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N 



(1) 
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where the spin operators are given in units of h, g is the 
spectroscopic splitting factor, and /is is the Bohr mag- 
neton. In the first term of Eq. the cyclic character 
of the system is maintained by requiring that sjv+i = sj- 
The second term describes the standard Zccman effect, 
where the external field H is typically defined to be di- 
rected along the z-axis. The total spin operators S 2 and 
S z then commute with H, and the eigenstates are de- 
scribed by quantum numbers S and Ms whose values 
range from to Ns and from —S to 5, respectively. In 
Fig. n we display the zero-field energy spectrum corre- 
sponding to Eq. for a particular example, s = 3/2 
and N = 6, with the subset of minimal energies (SME) 
shown in red. The SME are closely related to what are 
called level-crossing fields, quantities which are used to 




FIG. 1: (Color online) Complete energy spectrum for N = 6 
and s = 3/2 obtained by diagonalizing the Hamiltonain in 
Eq. Q for H — 0. The energy levels shown in red are the 
subset of minimal energies (SME) for this system as explained 
in the text. The solid line is a parabola, given in the text, 
that gives an excellent fit to the SME. 



study the SME in great detail throughout this work for 
many choices of s and N. 

In an external magnetic field, the (25 + l)-fold degen- 
eracy of each field-free multiplet is lifted due to a shift, 
gfisHMs, originating in the Zeeman term. As the ex- 
ternal field is increased from H = 0, the ground state 
will change (among the members of the zero-field SME) 
successively from 5 = 0, M s = to S = 1, M s = -1, 
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etc., in integer steps of S and Ms until S = —Ms = Ns, 
corresponding to saturation of the magnetization. Each 
of the Ns changes of the ground state quantum numbers 
is referred to as a level-crossing, and the field at which 
the nth level-crossing occurs is denoted in the following 
by H n . By determining these fields, we seek to record 
the characteristics of the SME as a function of s and N. 
This is accomplished using the difference equation, 

E N (S,s)=E N {S-l,s)+ gt i B H n (2) 

for S = n, where n extends from 1 to Ns. We elaborate 
on this connection between the SME and the H n in detail 
in the following section. 

In order to appreciate the details of the SME, we first 
review some generic features of the spectra, and in par- 
ticular the SME, that are already known. It has been 
noted ) 12 ' 1 ? and is clearly evident in Fig. ^ that the SME 
are accurately approximated by a quadratic dependence 
on S of the form E(S) cx S(S + 1), as for a quantum 
rotor. The solid curve in Fig. ^ describes the parabola 
E(S) = cJS{S + 1)/(2N) + E G , where c = 4.14 gives 
the best fit to this SME which has a ground state en- 
ergy Eq = — 17.393J. [The reason for the inclusion of 
the factor of 27V in this equation will become clear in 
the next section.] If the SME were strictly parabolic in 
S 1 , this would give rise to uniformly spaced level-crossing 
fields. Although uniform spacing is approximately real- 
ized in Fig.[5]for our example, we find that the accuracy 
of such an approximation deteriorates for larger values of 
N. This is explored in detail in Sec. ITD 

Above the SME there exists a large forest of energy 
levels. Although many of these levels lie very close to 
one another, there is a relatively large energy separation 
between the SME and the higher energy levels which has 
been previously observed^ii& Since at a low tempera- 
ture T only the lowest levels can be thermally occupied, 
and all other levels lie well above the SME, the mag- 
netization as a function of field M(H) consists of a sc- 
ries of thermally broadened steps that arise at the level- 
crossing fields and are determined solely by the SME. 
[The magnetization is also a function of T, but we will 
write M{H) for the sake of brevity] This step-like prop- 
erty is illustrated in Fig. El where M (if) and dM/dH, 
the differential susceptibility, are shown for the s = 3/2, 
N = 6 example with k B T / J = 0.1. The data in Fig.EJa) 
were calculated in three different ways: using the par- 
tition function that includes the complete energy spec- 
trum; using the quantum Monte Carlo (QMC) method 
employed in this work; and using a modified partition 
function that includes only the states belonging to the 
SME. The sharp peaks that appear in Fig. H[b) were 
calculated using QMC and the susceptibility fluctuation 
formula to give dM/dH directly, not by differentiating 
M{H). 

The three data sets shown in Fig.[3fa) are all identical 
to 4 significant figures, supporting the assertion that the 
SME are sufficient for analyzing low temperature exper- 
imental data of this type. For larger values of N, lower 




FIG. 2: (Color online) (a) M(H) and (b) dM/dH for the ex- 
ample system (N = 6 and s = 3/2) at a fixed temperature, 
ksT/J = 0.1. The data in (b) are obtained by QMC calcu- 
lations, and as a guide to the eye the line connects successive 
data points. 



temperatures are needed in order to obtain this degree 
of agreement, especially in the vicinity of the saturation 
field. For this reason, we have carefully checked that 
as the temperature is lowered the level-crossing fields 
have indeed converged to their limiting, temperature- 
independent values. As shown in Sec. |H] it is these 
fields that are then used to calculate the SME function 
E N (S,s). 

Despite the very simple appearance of the Hamilto- 
nian of Eq. Q , the evaluation of the corresponding en- 
ergy eigenvalues and resulting thermodynamic proper- 
ties frequently presents a major challenge. The most 
straightforward way to deal with this Hamiltonian, and 
the method that is usually employed when analyzing 
magnetic molecules, is to numerically diagonalize the 
Hamiltonian matrix. This yields energy spectra such as 
that shown in Fig. ^ However, even for relatively small 
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rings the dimensionality of the Hilbcrt space, given by 
D = (2s + 1) N , is so large that the exact diagonalization 
of the Hamiltonian matrix becomes totally impractical. 
For the small ring that has been considered as an exam- 
ple, s = 3/2 and N = 6 gives D = 4096. If we consider a 
larger ring, for example s = 5/2 and N = 12 which will 
be analyzed in Sec. IIVI we already have D ~ 2.2 x 10 9 , 
which is well beyond the practical limit of existing com- 
puters. For s = 5/2 and N = 20, D is a staggering 
3.6 x 10 15 . 

We can entirely avoid the obstacles confronting ma- 
trix diagonalization by using a QMC method that is not 
restricted by the dimensionality of the Hilbert space. 
We here only focus on low temperature M(H) and 
dM/dH(H) which are used to determine the SME, but 
other thermodynamic quantities such as the temperature 
dependent susceptibility, specific heat, and internal en- 
ergy are also readily attainable using this method for all 
temperatures and fields, and are in fact computationally 
much less demanding than the present low temperature 
studies. 

As seen above, knowledge of the SME enables one to 
obtain accurate values of low-temperature M{H) and 
dM/dH(H) data. To this end, the SME are calculated 
in Sec. [H] for all s from 1/2 to 5/2 and all even N from 
4 to 20. These data are presented in the form of con- 
venient, dimcnsionlcss "spectral coefficients" that will be 
introduced in that section. The spectral coefficients are 
also presented for larger rings, N = 40, 80, and larger 
intrinsic spins, s = 3 and s = 7/2. Such large values of 
N and s have not yet, to our knowledge, been realized 
in magnetic molecules, but are useful for studying the 
approach to the classical limit (s — > oo). 

In Sec. IIIII the energy gap A a (N) between the 5 = 
ground state and the lowest 5 = 1 state, which can be 
inferred from the first level-crossing field, is analyzed in 
greater detail for successively larger values of N, up to 
N = 50, for s = 3/2,2,5/2. This gap is experimen- 
tally relevant for NMR and INS experiments, and is also 
important for analyzing low temperature, low-field sus- 
ceptibility data. Finally, as an illustration of the use- 
fulness of the present results, in Sec. IIVI we analyze an 
existing ring-type magnetic moleculaii composed of 12 
Fe 3+ ions (s = 5/2), leading to an improved estimate 
for the exchange constant. With the experimental ad- 
vancements that are being made both in the synthesis of 
molecules and in high field magnetization studies, we an- 
ticipate that the use of the theoretical data presented in 
this work will complement future experiments in a much 
needed way, providing more accurate estimates of micro- 
scopic parameters for future ring- type molecules. 



II. SPECTRAL COEFFICIENTS 

Since the Hilbert space associated with H is often too 
large to allow diagonalization of the Hamiltonian matrix, 
other theoretical methods must be found. To analyze 



low-field susceptibility data, classical spin models and 
scaled-up data from smaller systems can sometimes be 
uge f u l ( 6 4 ii i i2 4 i4 jj oweverj the level-crossings that are ob- 
served in high-field experiments have no classical analog 
and cannot be easily scaled up. For this reason, reli- 
able theoretical data has previously been lacking, and a 
main goal of the current work is to remedy this situation 
through detailed QMC calculations. 

In order to learn about the nature of the SME, we used 
the Stochastic Series Expansion methodi2*22i to simulta- 
neously calculate both M and dM / dH versus H at fixed 
temperatures, an example of which was shown in Fig. [21 
From these data, we can very accurately infer the level- 
crossing fields and thereby reconstruct the SME. This 
follows from Eq. @ which gives 

s 

E N (S,s)=gn B J2Hn+E G , (1 < 5 < Ns), (3) 

n=l 

where Eg = En{0,s) is the ground state energy. It is 
convenient to define the quantities, 

_ gfj, B H„ c n {N,s)n 
h„ = — - — = — , {n=l,...,Ns), (4) 

where the dimensionless numbers c„ (N, s) will be referred 
to as "spectral coefficients" . The energy spectrum of the 
SME may thus be written as 

j s 

E N (S, s) = -Y, nc "( N > s ) + e g, (1 < S < Ns). 

(3') 

Note that if c n {N, s) were independent of n and given by 
c(N,s), Eq. Jsjij would reduce to 

c(N, s )JS(S + l) (n * q *„A 
E N {b,s) = — VE G , (0<b<Ns), 

(5) 

the so-called "rotational band" model that has often 
been employed to analyze magnetization dataAiLSLSS 
Inspecting Eq. the rotational band model imme- 
diately implies that the level-crossing fields are equally 
spaced which, as we will demonstrate in the subsequent 
subsections, is hardly the case. Instead, Eq. p 7 ]). in con- 
junction with the spectral coefficients presented in the 
following subsection, provides a highly accurate, yet con- 
venient means of representing En(S, s) and thus for an- 
alyzing low temperature magnetization data. 

Based on previously known properties of Heisenberg 
rings, it is easy to show that c n (N,s) is exactly 4 for 
a very few special cases. These are listed here and will 
be useful in discussing the results of our calculations in 
subsequent subsections: 

I. In the case of the N = 4 ring, c„(4, s) = 4, in- 
dependent of n and s. This is easily derived by 
describing this system in terms of two sublattices, 
each consisting of two spins. As a result, the SME 
is given exactly by E 4 (S, s) = J5(5 + l)/2 + E G . 
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II. In the limit of classical spins, 42 lim s ^ oc c n (N, s) = 
4, for all n and N. 

III. In all cases c n= N S (N 7 s) = 4, independent of N 
and s. This follows from Eq. ^} and the fact that 
the state with S = Ns has energy En(Ns,s) = 
JNs 2 , while the SME energy with S = Ns - 1 is 
E N (Ns -l,s) = Js(Ns - 4)M. 

Since the spectral coefficients have a value of exactly 
4 both in the limit of very small rings (item I) and in 
the limit of very large intrinsic spins (item II), one might 
expect that the replacement, c n (N 7 s) = c(N 7 s) ps 4, 
independent of n, would provide a very good approxi- 
mation, for example, for Fe 3+ ions (s = 5/2) in small 
rings (N < 20). However, as shown in the following 
subsections for different choices of s and N, the spectral 
coefficients do vary significantly with n. 



A. s = 1/2, and s = 1 

Rings of s = 1/2 spins have been studied using many 
methods, and a great deal is known about their spec- 
tra. In the 1960s, the lowest energies, En =00 (S,s = 
1/2), were calculated^ in the thermodynamic limit using 
the Bethe ansatzj2£i2£ while numerical diagonalization— 
was carried out on finite rings. More recently, work 
has continued for finite N using methods that include 
the quantum Monte Carlo2&, renormalization group2fli2i, 
Lanczosi&S and conformal field theory methodsi2i24 

The lowest eigenvalues for small s = 1/2 rings can be 
easily obtained from straightforward matrix diagonaliza- 
tion, but are included here both for completeness and to 
assess the usefulness of Eqs. I^J and (JHJ. The spectral 
coefficients that are shown in Fig. Efa) as a function of 
n/{Ns) define the SME for small s = 1/2 rings. One can 
immediately notice that the c n vary with n, and most 
are much larger than 4, implying that a rotational band 
approximation provides a relatively poor approximation 
to these spectra. 

Also included in Fig. [5{a) are the spectral coefficients 
corresponding to Griffith's original M(H) result for the 
infinite s = 1/2 chainS^ which is shown as a solid curve in 
Fig-HHb). In the thermodynamic limit, the transforma- 
tion from magnetization to spectral coefficients can be 
accomplished by making the replacement, n/N —> mo, 
where mo is the zero temperature magnetization per spin 



For rings of classical spins with N even, the SME can be de- 
scribed by the continuou s fu nction Epf(S c ) = 2J C S^/N + Eq, 
given in Eq. (80) of Ref. I24t Replacing the classical exchange 
constant J c = s(s + 1) J, and the classical spin S c by their quan- 
tum analogs, we obtain En(S,s — » oo) = 2JS(S + 1)/N + Eq 
from which item II follows. 
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FIG. 3: (Color online) (a) Spectral coefficients for small 
s = 1/2 rings. The solid lines are included to guide the eye, 
and the continuous curve corresponds to the N = oo magne- 
tization data of Ref.|2l (b) M(H) for N = oo (from Ref.|2j|, 
for N — 20 (QMC) and an approximation based on Q. 

in units of gfJ-B- Eq. Q can then be rewritten, 

c n (N = oo, s) = — ^— . (6) 
m {h n ) 

As can be seen in Fig. [3Ja) , for N = oo the spectral 
coefficients form a nearly linear function of n/(Ns) over 
a very large range of this variable. Approximating this 
data as a linear function, 

c n (N = co.1/2) w a -/?— , (7) 

and substituting Eq. (JJJ into Eq. again replacing 
n/N — > mo, the resulting approximate magnetization is 

where s = 1/2. Fitting the c„(oo,l/2) data to the lin- 
ear function, we find a = 8.9 and (3 — 5.07. The cor- 
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FIG. 4: (Color online) Spectral coefficients for small rings of 
intrinsic spins s = 1. The solid lines are included to guide the 
eye. 



responding curve terminates at the point (h = 1.953, 
mo = 0.439), rather than at (2, 0.5), but otherwise is 
virtually indistinguishable in Fig. EJb) from the exact 
magnetization (solid curve). This deviation of the termi- 
nus is due to the fact that the linear approximation of 
Eq. J7J) does not incorporate a small positive curvature of 
the cluster coefficients as a function of n/(Ns) as n ap- 
proaches Ns. Also included in Fig.^b) is M(H) for the 
s = 1/2, N = 20 ring at a temperature k B T/J = 0.05. 
This data is nearly identical to that of the infinite ring, 
except for the existence of thermally broadened steps as- 
sociated with level-crossings. 

Heisenberg rings of s = 1 spins have received a great 
deal of attention since Haldane's prediction^ that a finite 
gap separates the ground state from the first excited state 
in infinite rings of integer spins In the notation of the 
present work, this gap is given by A S (7V) = En(1,s) — 
Eg = jjCi(N, s), and the values of Ci(N, 1), seen as the 
left-most points in Fig. 0] are in good agreement with 
published values^! of A 1 (N). Values of A S (N) for all s 
in the range 1/2 < s < 5/2 will be discussed in Sec. IIIII 

Note however that the data presented here and in the 
next subsection include not only the first energy gap [as- 
sociated with Ci(N, s)], but all energy levels that be- 
long to the SME. Studying the details of of the SME, 
we find a very rich structure. For instance, it is evident 
in Fig. 0] that c n (N, 1) decreases rapidly with increas- 
ing n, unlike the corresponding data for ,s = 1/2. For 
?i > 0.4iVs the value of c n (N, 1) has already fallen be- 
low 5 for N < 20, whereas for s = 1/2 this value is not 
reached until n > 0.75Ns. In this sense, increasing s 
from 1/2 to 1 is a significant step on our way toward the 
classical limit, stated in item II of the previous section. 
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FIG. 5: (Color online) Spectral coefficients for small rings of 
intrinsic spins s — 3/2, 2, 5/2. The solid lines are included to 
guide the eye. 



3/2, 



2 and 



5/2 



Systems of larger intrinsic spins have also been studied 
in recent years jiLiLiSiiSiiSi but with less frequency than 
s = 1/2 and s = 1 systems. Since a knowledge of the 
spectral coefficients for s = 3/2, 2 and 5/2 is important 
for a number of molecular rings, these data are presented 
in Fig. for all TV < 20. The values of Ci(iV,s) that 
appear in Figs. EJa) and EJb) agree with the values of 
A S (N) that have been published^ (N < 10). Again, be- 
sides the first energy gap, the SME exhibit several inter- 
esting characteristics which arc reflected in the spectral 
coefficients. 

Of course, the spectral coefficients for N = 4 are all 
equal to 4 as required by item I. As N increases with fixed 



2 4 6 8 10 

gn B H/J 

FIG. 6: Differential susceptibility for a ring of N = 20 in- 
trinsic spins s — 5/2 at a temperature, ksT/J = 0.05. The 
large peak immediately before saturation (gfiBH/J — 10) is 
discussed in the text. 

s and n/(Ns), the corresponding spectral coefficients in- 
crease from 4 monotonically, resulting in the series of 
noninterseeting curves seen in Fig. [S] This is consistent 
with Waldmann's observation^ that the rotational band 
model becomes poorer for larger rings. 

Anchored at 4 for n = Ns (item III), and always ap- 
proaching 4 from above, the values of the spectral co- 
efficients decrease sharply as n approaches Ns. This 
ubiquitous drop can be discussed in a number of con- 
texts. Recalling Eq. |J2J, this is clearly equivalent to a 
compression of the level-crossing fields as saturation is 
approached. At low temperatures this results in a large 
slope of M(H), as can be seen in Fig. H3 for N — 20 and 
s = 5/2. In terms of the energy spectrum, this implies 
that the curvature of the SME decreases for large S. 

Finally, note that as s is increased with fixed N and 
n/(Ns), the spectral coefficients descend toward 4 (item 
II), but only very slowly. Even for s = 5/2, most of the 
spectral coefficients shown in Fig. OJc) are considerably 
larger than 4, indicating that one is still far from the 
classical limit that is stated in item II. This behavior 
is explored in the next subsection with the inclusion of 
larger values of intrinsic spin. 



C. Scaling relation for large s 

Thus far we have presented the spectral coefficients 
that define the SME as a function of three variables, s, 
N, and n/(Ns), and some general trends have emerged. 
Now, considering larger values of s and N, we would 
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FIG. 7: (Color online) Spectral coefficients adjusted to inves- 
tigate the scaling behavior of Eq. for all s in the range 
3/2 < s < 7/2. The data shown correspond to p = 1.05. 



like to make more quantitative statements regarding the 
functional dependence of c n (N, s) on these variables. To 
that end, we have calculated the spectral coefficients for 
values of s up to 7/2 and present that data for 3/2 < s < 
7/2. 

In Fig. H we plot the quantity [c n (N 7 s) — 4] x s p as a 
function of n/{Ns) for the choice p = 1.05. From these 
data the s dependence of the spectral coefficients is im- 
mediately evident. For each value of N, the data for all 
s lie on a single curve, implying that the spectral coeffi- 
cients scale according to 

c n (N, s) = 4 + s~ p f(N, n/(Ns)). (9) 

In particular, for s — > oo Eq. Q will be in accord with 
item II. The slow approach to 4 as s is increased is note- 
worthy, as even s = 7/2 is still far away from the classical 
limit. Choosing a slightly different value for the scaling 
exponent p, such as 1.03 or 1.07, yields visibly inferior 
results, so we conclude that p = 1.05 ± 0.01. 

A few of the spectral coefficients are also calculated 
for larger rings, N = 40 and N = 80. The inclusion 
of this data in Fig. [7| serves two purposes. First, this 
data suggests that f(N,n/(Ns)) is indeed converging to 
a finite limiting curve in the limit N — > oo, which defines 
the zero temperature M(H) of an infinite chain of spins 
s. Secondly, the larger N data strengthens our belief that 
the scaling relation © is valid for all N. 

Note that in Fig.|7|data are only included for n/ {Ns) > 
0.1. The data for small n/(Ns) have not been in- 
cluded because the error in calculating c n using the QMC 
method rapidly increases as n/(Ns) decreases towards 
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zero. The n 
section. 



1 (gap) behavior is considered in the next 



III. ENERGY GAP 

We now explore the energy gap A B (N) between the 
ground state and the first excited SME level. Values 
of this gap are shown in Fig. 03a) for rings of N < 20 
spins s < 5/2. Much like the behavior of the full SME 
discussed in the previous section, this gap systematically 
approaches the limiting s — oo form as s increases from 
3/2, while the s = 1/2 and s = 1 data exhibit distinctly 
different trends. 

Specifically, the energy gap for s = 1/2 rings is very 
similar to the energy gap that would be obtained for a 
ring with the same value of N but very large s. This 
large s limit, indicated in Fig. EJa) as the "classical ro- 
tational band" , follows from item II and is given by 
Aoo(TV) = AJ/N. By contrast, s = 1 rings have much 
larger gaps. Note also that these are already within 3.5% 
of the limiting, N = oo value even for N = 20. The 
known limiting valued Ai(oo) ~ 0.4105 J, is indicated 
by an arrow on the right side of Fig. 0H a )- 

Recall that for any s, A s (4) = J from item I. Consid- 
ering N = 6, the classical result Aoo(6) = 2 J/3 is still a 
reasonable approximation to A s (6), with a relative error 
of only a few percent for any s. However, as N increases 
further this error continues to grow, and with N = 20 it 
is nearly 25% for s = 5/2 and nearly 40% for s = 3/2. 

Although the classical result is not sufficient, we find 
that the energy gaps for s > 1 are well described by a 
slightly more general power law dependence on N of the 
form 



A S (N) ~ riN~ 



(10) 



The curves in Fig. Ufa) were obtained by choosing: fl = 
2.68 J and a = 0.757 for s = 3/2; Q = 2.73 J and a = 
0.781 for s = 2; fl = 3.03 J and a = 0.837 for s = 5/2; 
while of course 51 = 4 J and a = 1 corresponds to s = oo. 
With these choices of f2 and a, excellent agreement with 
the QMC data is obtained in the range 8 < N < 20, 
and it is clear that the classical limit is indeed being 
approached with increasing s for both fl and a. 

For the half odd integer spins, s = 3/2 and 5/2, the 
agreement with Eq. (|10f> continues for larger values of 
N. The same data are shown in Fig. 02b), now includ- 
ing N < 50, and the QMC data agree with the power 
law formulas to within a fraction of a percent for all ring 
sizes in the range 8 < N < 50, which is comparable to 
our uncertainties in A S (N). The values of A2(N) begin 
to diverge from the power law dependence for N > 30 
which is to be expected since they must eventually con- 
verge to a non-zero value. This gap for s = 2 chains 
has been previously studied in great detail, and density 
matrix renormalization group calculations have yielded 
a valueii of A 2 (N) = 0.0876 J ± 0.0013 J in the limit as 
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FIG. 8: (Color online) Energy gap A B (N) for rings of N 
spins (a) for all s < 5/2 (symbols) along with the best fits 
to Eq. (llOt (curves) as described in the text; and (b) with N 
varying from 4 to 50 for s — 3/2, 2, 5/2 (symbols) along with 
the best fits to Eq. 1101 (lines) as described in the text. 



N — > oo. One can see in Fig.OJb) that A 2 (N) is begin- 
ning to approach its limiting value, having become larger 
than A 3 / 2 (A'') for N > 50, but data for much larger rings 
would be necessary in order to obtain an accurate esti- 
mate for the limit N — > oo. 

The rotational band result, A S (N) = AJ/N, has been 
used in the pastiii*i2i2i as an estimate of A S (N). Al- 
though this provides a reasonable approximation for N < 
10, as we have seen it quickly diverges from the correct 
result with increasing N. As such, it would be prudent 
to use the more accurate results presented here when 
attempting to relate J to the experimentally measured 
lowest energy gap, e.g., by using INS, NMR, low-field 
susceptibility or magnetization data. 
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FIG. 9: (Color online) The four measured level-crossing 
fields— (dashed vertical lines) are compared with the theo- 
retical dM/dH that result from N = 12, s = 5/2 Heisenberg 
rings with ksT/J = 0.01. The theoretical data are shown 
for the following three choices of J and g: J — 31.9 K and 
g = 2.00 (o) from Ref. [2 J = 40.7 K and g — 2.02 (□) from 
Ref. l2lt J = 35.2 K and g — 2.0 (o) are our best estimates. 



IV. AN APPLICATION: Fei 2 

In this section wc apply our results to a known mag- 
netic moleculejii whose analysis has been challenged by 
a Hilbert space of dimension D = 6 12 ps 2.2 x 10 9 . The 
molecule is comprised of 12 Fe 3+ ions (s = 5/2), whose 
interaction was first investigated^ by measuring the low- 
field susceptibility Xo{T) as a function of temperature 
and fitting this data to an approximation of the s = 5/2, 
N = oo chain. The exchange energy thus obtained was 
J/k B = 31.9 K for g = 2.00. The field dependent mag- 
netization of the molecule has also been measured and 
analyzed, and the first four level-crossing fields at low 
temperatures wereSI Hi = 10.1±0.2T, H 2 = 19.6±0.2T, 
H 3 = 29.6 ± 0.4T, H 4 = 39.1 ± 0.8T. An analysis of the 
magnetization was given in Ref. |^] using the classical ro- 
tational band c(N, s) = 4, and this yielded the estimate 
J/k B = 40.7 K with g = 2.02. Note that the latter esti- 
mate is more than 25% larger than the former— estimate. 
Given the results of Sec. [H] one can expect that the es- 
timate J = 40.7 K will be considerably larger than what 
will result from an accurate treatment of the Heisenberg 
model. This is borne out in the following. 

In Fig. [5] wc compare the four measured level-crossing 
fields with our low temperature (k B T/J = 0.01) QMC re- 
sults. At this low temperature, each level-crossing of the 
theoretical s = 5/2,A = 12 Heisenberg ring is clearly in- 
dicated by a narrow peak in dM/dH . Note that the peaks 
in the QMC data arising from the parameters J = 31.9 K 
and g = 2.00 occur at fields that are considerably below 
the experimental level-crossings indicated by the dashed 




FIG. 10: The first four spectral coefficients for the N = 12, 
s = 5/2 ring (x) along with their errors bars, compared 
with the results of inserting the first four level-crossing fields 
and uncertainties from Ref. |^ into Eq. (|1J with the choice, 
J/(k B g) = 17.6 K (o). 



vertical lines. On the other hand, the QMC peaks that 
correspond to J — 40.7 K and g = 2.02 are all at fields 
significantly greater than the measured values. Particu- 
larly at high fields, these discrepancies become quite pro- 
nounced, suggesting that neither choice of parameters is 
consistent with the experimental data. However, we find 
that the predictions of the Heisenberg model agree very 
well with the experimental data if we use J = 35.2 K 
and g = 2.0. With this choice of parameters, each of the 
four theoretical peaks clearly coincides with a measured 
level-crossing shown in Fig. [5] 

Without using the dM/dH level-crossing field data di- 
rectly, one can easily arrive at the same estimate based 
on the spectral coefficients of Sec. |H] Recalling Eq. (0J, 
the ratio of J to g is given by J/g = Nfj,BH n /[nc n (N, s)]. 
An estimate of this ratio for a given molecule is then ob- 
tained by simply inserting a measured value of H n and 
the corresponding c n (N,s) from Fig.|SJ 

Alternatively, from the measured H n we can construct 
an experimental analog of the spectral coefficients by fix- 
ing the ratio J/g in Eq. In Fig. ^2|wc display those 
results for the four measured H n (and their uncertain- 
ties). These data are in good agreement with the spec- 
tral coefficients if we choose the ratio J/(ksg) = 17.6 K, 
consistent with our previously stated estimate. 

A small decrease with increasing n is observable in the 
spectral coefficients derived from the experimental values 
of the level-crossing fields. This is expected from the 
data presented in Scc.[Hl but more level-crossings and/or 
smaller experimental error bars are needed in order to 
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clarify this point. These data are also useful for getting 
a sense of the typical errors in the spectral coefficients 
that were presented in Sec. |n] As shown in Fig. EH the 
error bars of the QMC data decrease very rapidly with 
increasing n and are in fact not visible in Figs. [5] and 
Our conclusion is that the existing data for the Fei2 
molecule is best fit by the choice g = 2.0, J = 35.2 
K. This value of J is 13.5% smaller than the value that 
resulted^ from analyzing the experimental level-crossing 
fields using c(N,s) = 4. This reflects the fact that the 
spectral coefficients, although not constant, exceed 4 by 
approximately 13%. A similar analysis would be equally 
straightforward for any other rings whose spectral coeffi- 
cients are shown in Sec. [H] 

V. SUMMARY 

In this article we have utilized a quantum Monte Carlo 
(QMC) methodiSiSS to calculate detailed properties of 
the general quantum Heisenberg ring. This system con- 
sists of an even number A of equally-spaced spins s 
mounted on a ring, where the spins interact via nearest- 
neighbor antiferromagnetic isotropic exchange, with a 
single exchange constant J. As this system does not 
exhibit magnetic frustration it was possible to calcu- 
late thermodynamic quantities down to very low tem- 
peratures. In this work our primary focus has been on 
the accurate determination of the level-crossing fields, 
which in turn directly provide the lowest energy eigen- 
value En(S, s) for each value of the total spin quantum 
number S. By introducing the notation of spectral coeffi- 
cients [see Eq. Q], denoted by c n (N,s), we obtained an 
especially convenient representation of En(S,s), given 
by Eq. J3]). As the QMC method operates without diag- 
onalizing the Hamiltonian matrix, we were able to obtain 
results for spins s = 1/2, 1, 3/2, 2, 5/2, 3, 7/2, focusing 
mostly on A < 20 as these are experimentally relevant, 
although N < 80 were also considered. Among our prin- 
cipal results, we have found that the set of level-crossing 
fields are not uniformly spaced, and thus the spectral co- 
efficients c n (N, s) vary significantly with n. Equivalcntly, 
Ejy(S, s) departs from the strictly quadratic dependence 



on S*, referred to as the rotational band approximation^ 
(equivalently, the Lande interval rule). These deviations 
from uniform spacing are fairly small for N < 8, however, 
they become increasingly severe with increasing A. Sim- 
ilarly, the ground state energy gap, which may be written 
as A S (N) = E N (l,s) - E N (0,s) = a(N,s)J/N, varies 
significantly with A and s. For s = 1, 2 we find that, 
consistent with the Haldane result^ A S (A) is indeed 
converging to a non-zero limiting gap for large A, and 
in good agreement with estimates in the literatures*^ 
for these two choices of s. By contrast, for s = 3/2, 5/2 
we find that A s (A) appears to decrease to zero for large 
A according to a power law, [see Eq. I|10[) ]. where the 
exponent is given by a k 0.76 for s = 3/2 and a « 0.84 
for s = 5/2. The increase with a towards unity with in- 
creasing s is consistent with the known rigorous result^ 
Aoo(A) = 4J/A, for the classical Heisenberg ring, which 
may be pictured as the quantum Heisenberg ring in the 
limit s — > oo. We also find that the departure of the 
general spectral coefficient c n (A, s) from the classical re- 
sult c n (A, oo) = 4 is characterized by power law behav- 
ior [see Eq. (JUJ)]. Finally, we have illustrated the prac- 
tical utility of our present results for the level-crossing 
fields and En(S,s) by considering the ring- type mag- 
netic molecule knownii as Fei2- In particular, our anal- 
ysis of the existingSi experimental data for level-crossing 
fields shows that this system can be very well described 
by the nearest-neighbor Heisenberg model with antifer- 
romagnetic exchange, and we are able to provide a new 
and improved estimate of the exchange constant. Al- 
though only rings with even A have been considered in 
this work, we suspect that similar scaling relations may 
hold for other structures as well. 
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